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Abstract 

Using a TE/TM decomposition for an angular plane- wave spectrum of free random elec- 
tromagnetic waves and matched boundary conditions, we derive the probability density func- 
tion for the energy density of the vector electric field in the presence of a semi-infinite isotropic 
medium. The theoretical analysis is illustrated with calculations and results for good elec- 
tric conductors and for a lossless dielectric half-space. The influence of the permittivity 
and conductivity on the intensity, random polarization, statistical distribution and standard 
deviation of the field is investigated, both for incident plus reflected fields and for refracted 
fields. External refraction is found to result in compression of the fluctuations of the random 
field. 

1 Introduction 

Complex electromagnetic (EM) environments are characterized by the fact EM fields behave 
as random or quasi-random quantities. They can be studied in an efficient manner with the 
aid of statistical electromagnetics methods. With regard to the interaction of random fields 



with their environment, a fundamental problem of interest is the evolution of statistical prop- 
erties of the field upon propagation through stratified media, including reflection and refraction 
at their interfaces. For a perfectly electrically conducting (PEC) surface, previous studies of 
the average value [Dunn, 1990], standard deviation and first-order probability density function 
(pdf) [Arnaut and West, 2006] of the electric and magnetic energy densities have demonstrated 
a direction-dependent damped oscillatory behaviour of their average value and standard devi- 
ation as a function of the distance of the point of evaluation to the interface. This behaviour 
is a consequence of the interference between incident and reflected fields. As a result, unlike 
for deterministic waves, a boundary zone exists for random fields adjacent to the PEC surface, 
in which the statistical field properties are inhomogeneous and fundamentally different from 
those at larger (theoretically infinite) distance. Further insights that were gained from these 
studies pertain to the statistical anisotropy and polarization state of the field within the bound- 
ary zone and, for vector fields, the transitions of the pdf of the energy density from one- or 
two-dimensionally confined random fields at the interface to fully developed three-dimensional 
random fields at large electrical distances. In addition, spatial correlation functions have been 
obtained previously for unbounded [Bourret, 1960], [Sarfatt, 1963], [Eckhardt et al], [Mehta and 
Wolf, 1964], [Hill and Ladbury, 2002] and single-interface [Arnaut, 2006a], [Arnaut 2006b] vec- 
tor EM fields that elucidate the spatial structure of random fields via their two-point coherence 
properties. 

In the present paper, the methods and results for statistical properties of random fields near 
a PEC surface are extended to a magneto-dielectric isotropic semi-infinite medium. Having an- 
alyzed the second-order spatial coherence and correlation properties for an impedance boundary 
in [Arnaut, 2006b], here we are again concerned with first-order statistical, i.e., local distribu- 
tional properties only. Based on previous results for nonlocal spatial coherencies of the electric 
field (E a (ri)Ep(r2)) (a, (3 = x,y,z), the polarization coefficient and pdf for the local energy den- 
sity are determined. Because of the single interface and isotropy of the medium, the polarization 
coefficient is degenerate, whence the pdfs are one-parameter compound exponential (CE-1) dis- 
tributions [Arnaut, 2002], [Arnaut and West, 2006]. However, unlike for a PEC medium, the 
angular spectra of reflected and refracted random fields exhibit directivity because reflection and 
transmission coefficients of plane waves for a magneto-dielectric semi-infinite medium depend 
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on the wave polarization and angle of incidence. We shall confine the analysis and results to 
the electric field E; corresponding results for the magnetic field follow without difficulty. Since 
we express results in terms of Fresnel reflection and transmission coefficients for an isotropic 
half-space, analogous results for multilayer strata are easily obtained from the listed integral 
expressions, on substituting with the appropriate coefficients. In particular, probability distri- 
butions of reflected and transmitted fields on either side of a single layer of finite thickness are 
easily computed. 

2 Field coherencies 

Since random EM fields are spatially and temporally incoherent, their fundamental statistical 
quantity is the energy density - rather than the field itself or its magnitude - which is obtained 
from the average field intensities (l-E^r)! 2 ) of the Cartesian complex field components E a . These 
intensities represent self-coherencies, which are special cases of mutual coherencies [Mehta and 
Wolf, 1964]. The latter were derived in [Arnaut, 2006b] for the present configuration. Here, the 
fact that no separate transverse locations need to be considered in order to obtain first-order 
statistics considerably simplifies the calculations. 

2.1 Incident plus reflected random field 

We consider a semi-infinite isotropic medium with permittivity e and permeability /i occupying 
the half-space z < (Fig. [I]). The incident random field at r in the region characterized by 
z > is represented by a statistical ensemble (random angular spectrum) of time-harmonic 
plane waves [Whittaker, 1902], [Booker and Clemmow, 1950], [Hill, 1998]: 



in which an exp(jut) time dependence has been assumed and suppressed. This incident field 
E l is assumed to be ideal random, i.e., any three complex Cartesian components, in particular 
(E l x , Ey, E l z ), are mutually independent and exhibit identical circular centered Gauss normal 
distributions. The direction of incidence for each plane-wave component (£ l ,H l ,k l ) is arbitrary 
within the upper hemisphere (Oo = 2tt sr) and is specified by azimuthal (0o) and elevational (#o) 
angles. Since the medium is deterministic, the incident and reflected fields for each individual 




(1) 
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plane wave are mutually coherent, despite being individually random. Hence, their recombina- 
tion in the region z > is governed by superposition of fields, rather than energy densities. The 
boundary condition does not affect the correlation between field components on either side of 
the boundary, because £ l is itself random. As a result, the field components remain mutually 
independent in the vicinity of the boundary. Since the Cartesian components of the incident 
field are circular Gaussians and because the medium is linear, the incident plus reflected vector 
field and its components are also circular Gauss normal with zero mean, but now with different 
standard deviations from those of the incident field owing to the boundary conditions. Thus, 
the three Cartesian components no longer exhibit identical parameters for their distributions. 

Following [Dunn, 1990] and [Arnaut, 2006a], we perform a TE/TM decomposition for each 
plane wave component of the angular spectrum with respect to its associated random plane of 
incidence ok l z, i.e., cj) = cpQ. As is well known, TE and TM polarizations constitute uncoupled 
eigenpolarizations for stratified media, hence the polarization of the outgoing wave is completely 
determined by that of the the incoming wave. As a result, the TE and TM contributions to 
the overall plane-wave spectrum can be calculated independently. Ensemble averaging then 
yields the TE and TM energy contents of the random field. Specifically, for the TE components 
incident at an angle 0o an d with electric field 

£ j_ exp ^— jk* • r) = £q cos i/j exp (J k$x sin 0q) 

x exp (jk z cos O ) ly, ( 2 ) 

the incident plus reflected electric field at r is 

£± exp(-jk • r) = £ exp (-jk g sin O ) 

x {[1 + r_i_(0 o )] cos (fco^COS 00 ) 

+i[l-r ± (0 o )]sin(A;ozcos0o)}, (3) 

with q = x cos <po + y sin <pQ, and where r_|_(0o) represents the Fresnel reflection coefficient for TE 
waves for a semi-infinite isotropic medium. Upon substitution of ([3|) into an expression similar 
to (P) for Ej_, followed by unfolding and integration, the associated average field intensity is 

i\E y (r)\ 2 ) = I yl + I y2 + I y3 + IyA, (4) 
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in which (•) denotes ensemble averaging of the plane- wave spectrum, with 

I y i = 2C f 1 \l -T ± (u)\ 2 sin 2 {k zu)du, 
Jo 

2C [ 1 1 + r ± (it) | 2 cos 2 ( k zu) du, 
Jo 

I* yA = j2C £[1-T ± (u)} [i + ri(u)] 
x sin (kozu) cos (kozu) du, 



ly2 



iy3 



where u = cos 00, C = (\£o\ 2 )/4, and 



T ± (u) 



rjku — r] J k 2 — k 2 + k 2 u 2 



rjku + rjo^J k 2 — k 2 . + k 2 u 2 
Similarly, for the TM components with incident electric field 



fjjexpf-jV-r 



—So sin ip exp (jkox sin 9q) 
x exp (jkozcos 9q) 

X [cOS^ola; + sin#olz] , 



we obtain 



(\E x (r)\ 2 ) = J2l xe , (\E z (r)\ 2 ) = ^I ze , 



with 



Ixl 


= 2C / 




Jo 


1x2 


= 2C j 




Jo 


1x3 


T* — 
— 1 x4 — 



1-TJu) 



u 2 sin 2 (kozu) du, 



u cos (kozu) du, 



J2C f 1 
Jo 



!- r l|( 



u 



l + TUu) 



hi 

Iz2 
Iz3 



xu 2 sin (kozu) cos (kozu) du, 

2C ^ 1 + V\\ (u) 2 (1 - u 2 ) sin 2 (k zu) du, 
2C ^ 1 - V\\ (u) 2 (1 - u 2 ) cos 2 (jfeozu) du, 

4* 4 = i2C £ [i + r,|(«)] [i-q(n) 

x (1 — u 2 ) sin (fcoztt) cos (kozu) du, 



where 



r/Jk 2 -kl + klv? -rjoku 

T\\(u) = v = . (17) 

7] J k 2 — kg + k\u 2 + rjoku 



For the tangential and overall (i.e., incident plus reflected) vector field, 



(I^W| 2 )= E E^ ( 18 ) 

a=x,y £=i 

and 

(l^(r)| 2 )= E E^, (19) 

a=x,y,ze=l 

respectively. For a medium for which \/j, r /e r \ — > 0, only the terms I a i remain nonzero. Con- 
versely, for |e r //x r | — > 0, only the terms 7 a 2 survive. 

2.2 Refracted random field 

For the field refracted across the boundary, 



(\E y \ 2 ) = 2C f 1 \T ± (u)\ 2 du, 
Jo 



(20) 



du, (21) 



(l^l 2 ) = 2C j\ T]l (u)\ 2 ^ (l - u 2 ) du, (22) 



with the TE and TM transmission coefficients given by 



T ± (u) = (23) 

77 /cu + 770 a/ fe 2 — fcg + k^u 2 



and 



W = / = = ■ ( 24 ) 

rjoku + rjJ k 2 — k$ + fcgn 2 

respectively. Thus, unlike for the incident plus reflected field, the intensity of the refracted field 

is homogeneous due to the absence of interference in the region z > 0. 



3 Energy density distribution 

For a single-interface configuration, the pdf of the electric energy density (S e ) = e( )(|-E| 2 )/2 
can be calculated based on knowledge of the polarization coefficient P e ^ for the statistically 



uniaxial electric field [Arnaut and West, 2006]. For the vector tangential field, the incoherent 
superposition (\E t | 2 ) = (\E X \ 2 ) + (\E y \ 2 } holds, because £ x and £ y belong to mutually orthogonal 
modes. Hence, 

1 3(\E z (k z)\ 2 ) 

Pe,13-Pe,23 - 3 " 2{\E(k z)\*) ' (25) 

in which E z = E3, and 

<|£(M| 2 > = ]T <|P Q (M| 2 > (26) 

= 2C £ [\T ± (u)\ 2 + m(u)\ 2 ] du. (27) 
The associated CE-1 pdf of S e follows as [Arnaut, 2002], [Arnaut and West, 2006] 

fs e ( s e) = 7is e exp(-ais e ) + 72 exp(-a 2 s e ) +7 3 exp (-a 3 s e ) , (28) 

where 

ai = a2 = (TT^am) = J^TisZ) ■ (29) 

03 = (i-2p! 13 )(5 e ) = (i)' (3o) 

3 

7i : 
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P e ,13(l+Pe,13)(S e ) 2 

4 

«5e)-3(5 c ,))((5 B )-(5 c ,))' 
2 - P e 2 13 + ^e,13 ~ 1 

3P 2 13 (1 + P e , 13 )(5 e ) 



(31) 



73 



«S e )-3(S e J) 
4P 2 3 - 4P e , 13 + 1 



2 ' 



(32) 



3P 2 13 (l-2P e , 13 )(S e ) 

-72- (33) 



MS ez ) 



((S e ) - 3(S ez )Y 

4 Special cases 

4.1 Good electric conductors 

For nonmagnetic good electric conductors, i.e., a 3> weo, e = eo, fJ- = fJ-o, it follows that 



77/770 ~ y / uje /(2a)(l±j), T±(u) ~ (r]u - r] ) / (w + ij ) andrii(u) ~ (v ~ Vo u ) / (iJ + Vo u ) because 



7 



6 ~ irrespective of 8q. As a result, ([S])-© become 



«»2 



«»3 



4C 



4C 



1 1 — cos(2kozu) 



du, 



2 Z" 1 1 + cos(2/c .zu) n , 
' — u du. 



1 + 

VO 



i; 4 = j4c^- f 

y T)0 JO 



r/* f l sm(2k zu) 



udu. 



(34) 
(35) 
(36) 



Throughout the range of integration, |1 + (rj/rio)u\ 2 = [1 + u^uje^/ (2a)] 2 + u 2 [ujeo/(2a)] ~ 1. 
Therefore, to good approximation, we have 



(I^WI 2 ) 



AC 



-J- 



m 
H 

Vo 
rj — rj 



2 2cos(2/c z) 
(2k z)i 

2 sinc(2fco-z) 



+ 



sinc(2A;o^) 
sine (2koz) 



(2k z) 2 
sinc(2kQz) cos(2kQz) 
(2k z) 2 {2k z) 



(37) 



For the TM components, — (jl6|) specialize to 



'.xl 



1x2 



AC 



AC 



1 1 — cos(2kozu) 4 



u + 



'/() 



u du, 



2 rl 1 + cos(2k zu) n , 
— u du, 



u + 



'/() 



I x3 = I* x4 = jAC^- f 
rjo Jo 



rj* f 1 sm(2kozu) .. 
~ u 

u + ^ ' 



u 3 du, 



,^ f 1 1 + cos(2k zu) / 2 \ 2 , 
I z \ = AC / — 1 - u ) u 2 du, 

VO 

2 nl 1 — cos(2kozu) / 2 \ , 
— II — u J du, 



h-2 



'z3 



AC 



I 



zA 



jAC 



Vo Jo 



U+ V0 

sin(2kQZu) 



u + 



(l — u 2 ^j udu. 



(38) 
(39) 
(40) 

(41) 
(42) 
(43) 



Before presenting results for the pdfs, we show in Fig. [2] the dependencies of (S et ) and (S £z ) 
on koz at selected values of a/(ueo), after normalization by the average electric energy density 
of the incident total (vector) field, (S eo ). The sensitivity to variations of a/(ueo) is seen to 
be significantly higher for the normal component than for the tangential component. For finite 
a/(u€o), the value of (S et (koz— >0)}/{S eo ) is close to, but different from zero on account of the EM 



boundary condition. Both asymptotic values (S et (koz^+oo)} / (S eo ) and (S ez (fco-z— ►+oo))/(S eo ) 
are smaller for a finitely conducting boundary than for a PEC boundary, on account of energy 
dissipation in the former. For intermediate values of koz, finite conductivity also gives rise to a 
relatively small positive or negative phase shift in the oscillatory behaviour of (S ez (koz)) / {S eo ) 



The effect of finite conductivity on the statistical polarization state of the field is shown in Fig. 
[3l Unlike for a PEC surface, P^k^z) = P e ,i3(koz) (i = 1, 2) is no longer oscillating symmetrically 
with respect to zero, i.e., the random polarization exhibits a conductivity-dependent bias in the 
normal direction. This bias decreases, on average, with increasing koz, but persists nevertheless 
up to asymptotically large distances k$z — ► +00. 

Figure H] shows scaled pdfs fsX s e) for selected values of ojioje^) at koz = tt/4. At this 
intermediate distance, the pdfs rapidly approach the asymptotic distribution for S e near a PEC 
wall [Arnaut and West, 2006] when a/(ojeo) is increased. For larger k$z (not shown), the 
influence of a/(coeo) on the pdf was found to be even weaker. On the other hand, very close 
to the boundary {k^z <C 1), the pdf rapidly approaches a distribution close to the x| limit 
distribution when g/{uj€q) is increased, as demonstrated in Fig. [5] for k$z = 0.01. 

As indicated by Fig. [6|a) , the standard deviation (std) as e for the incident plus reflected total 
field exhibits increasing oscillations for intermediate distances and lower values upon approaching 
the surface when conductivity is decreased. The std approaches a location-dependent asymptotic 
value in the manner shown in Fig. [f^b). 

4.2 Lossless dielectric medium 

For a lossless medium, 77, k, T^y and F^y are real- valued. Therefore, and I a 4 are purely 
imaginary so that only I a \ and I a 2 contribute to the intensities (\E a \ 2 ). The following explicit 
expressions are obtained for the energy density of the fields refracted by a lossless dielectric 
medium with relative permittivity e r = e/eg: for TE wave components, 



or (S et {k z))/(S eo }. 




2C 



3(e r -l) 15(e r -l) 




(44) 
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whereas for the TM components, the integrals 



with 



and 



with 



Ui 

U 5 



V 2 
V 3 
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f 




du 


10 







2CJ2U 



4e r (e 2 - e r + 1) - 4e r 3 / 2 (2e r - 1) 

( er -l)2( £r + l)2 

12e 2 



(^-1)3/2^ + 1)2' 

2e r (2e 2 + 1) 



In 



(e r - l) 2 (6 r + 1)V2 

2e r (2e 2 + l) 
- l) 2 (e r + 1)V2 

2e r (2e 2 + 1) 
(e r -l)2(e r + l)5/2"M er + ^TZTJ 



In 



In 



v / e7TT- 1 \ 
VeTTT+lJ ' 

y/e r + 1 + 

VeTTT- y /e r ' / 
e r 



1 







2 






du 


10 





2C^V-, 

i=l 



4e r (e 2 + 1) 



8ef 



3/2 



(er - 1) 



3/2 



3( £ 2 - 1)5 



4e r (3e 2 - e r + l) - 4er /2 (4e r - 1) 
(er-l) 2 (e r + l) 3 
20e 2 



(e r - if 2 (€r + 1)' 
2e r (4e 2 + l) 



In 



(e r - l) 2 (e r + l) 7 / 2 ^ ^ v / e7TT + 1 
2e r (4e 2 - e r + l) 



( er _l)2 (er + 1) 7/: 



-In 



2e r (4e 



e r + l 



rln 



vV + 1 + ^/e7 

: er-Vg"£l \ 
e r + Ve 2 - 1/ 



(45) 



(46) 
(47) 
(48) 
(49) 
(50) 



(51) 



(52) 

(53) 
(54) 

(55) 

(56) 

(57) 



( er _l)2 (er + 1 )7/2 

can be substituted into (|2ip and (|22p to yield (l-E^I 2 ) and (lE 1 ^ 2 ) for the refracted field. 

Figure [7] shows the intensities of the tangential and normal components as a function of 
permittivity, for the incident plus reflected field at koz = 20n, and for the refracted field at 
any k$z. For the former case [Fig. da)], the minimum in the curves is a result of reduced TM 
contributions near the Brewster angle for external refraction (e/eo > 1). This effect occurs in 
both the tangential and normal field components, via (\E X \ 2 ) and (\E Z \ 2 ). However, for (\Et\ 2 ), 
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it causes only a minor dip in the characteristic because of the dominance of the contribution 
by (\E y \ 2 ) over that by (\E X \ 2 ). Note the large sensitivity to the permittivity in the vicinity of 
e/eo = 1. This sensitivity decreases with increasing Uqz. For the normal field, the minimum 
intensity occurs around e/eo — 5.4. The incident energy density for this component is being 
exceeded only above relatively high permittivities (e/eo ~ 70). For the refracted field [Fig. 
Efb)], the intensities rapidly decrease with increasing e/eo an d exhibit an increasingly dominant 
tangential contribution. 

Figure [8] shows the polarization coefficient as a function of e/eo- For the incident plus reflected 
field, an e r -dependent threshold distance exists where starts to increase, in an oscillatory 
manner, toward an (e r -dependent) asymptotic value when k$z — ► +oo. 

Figure [9] shows fs e (s e ) for the refracted field at selected values of e/eo- This pdf is indepen- 
dent of koz but evolves from a x| distribution for e/eo = 1 to a x\ distribution for e — ► +oo. 
From Fig. [TUJ <rs e (e/eo) is seen to increase monotonically from its x| value, l/\/3, at e/e = 1 to 
its asymptotic %\ value, l/v2, represented by the dotted line. 

Figure [TT1 shows fs e (s e ) for the incident plus reflected field close to the surface (koz = 7r/128), 
at selected values of e/eo- The pdf is seen to make an excursion as e/eo increases, returning 
eventually to the asymptotic x\ pdf when e — ► +oo. For larger values of k^z, it has been found 
that the excursions are much shorter; for example, for k§z ~ 2ir, there is essentially no longer a 
discernable parametric dependence of fs e (s e ) on e r . 

5 Effect of conductivity and permittivity on the magnitude of 
field fluctuations 

Since the probability distribution of each Cartesian component of the underlying complex field 
E a = E' a — jE'a remains circular Gauss normal, on either side of the interface with an isotropic 
medium, the standard deviation of E a can be deduced from the \\ statistics for | E a \ 2 or S Ca 
via 



of (JE a on a/(u>eo), e/eo and koz therefore follows immediately from the results shown in Figs. 




(58) 



where a = x,y or z, with (S ex ) = (S e ) = (S et )/2 and (je x = oe v = &E t /V2- The dependence 



11 



[2] and Physically, after re-scaling, Fig. E^b) indicates that an increased permittivity causes 
the fluctuations of the refracted fields to exhibit a smaller spread compared to the fluctuations 
of the incident field, where the latter are quantified by o~E a with e/eo = 1. This compression of 
fluctuations is seen to be more prominent for the normal field component than for the tangential 
field. By contrast, Fig. 0(a) indicates that the spread of the incident plus reflected field increases, 
on average, with increasing permittivity. 

6 Applications 

Besides its interest as the solution to a fundamental problem in statistical electromagnetics, the 
above analysis is relevant to several practical applications, of which we give a examples. 

• Atmospheric and ionospheric propagation of EM waves with applications, including spec- 
troscopy of stellar light: Stellar light propagating through an atmosphere can be repre- 
sented as a random field radiated by a collection of incoherent point sources that are 
distributed across one or several narrow solid angle(s) (cf. [Arnaut, 2006b] for a treat- 
ment of random fields produced by a spatially filtered EM beam). Upon refraction by 
atmospheric layers exhibiting permittivities close to eo, the light undergoes a change in 
its distributional and statistical polarization properties; cf. Figs. E^b) and [8jb). Fur- 
thermore, recall that the abovementioned plane-wave expansion refers to a harmonic or 
quasi-harmonic (narrowband) field for a specific central wavelength. Therefore, the pdfs of 
light intensities associated with different spectral lines may undergo different changes upon 
propagation, as a result of frequency dispersion of the medium, thus producing nonuni- 
form changes of the mean values and standard deviations across the received spectrogram. 
Similar considerations apply to scattered or reflected fields. In this way, the presented 
analysis yields corrections that could be instrumental in astrophysical observations and 
Earth sciences, including remote sensing. 

• Reflection and refraction of multipath signals for wireless communications: An important 
issue in the accurate prediction of multipath and wideband propagation through a radio 
channel is the determination of transfer functions for incident and outgoing signals, and 
how this transfer might affect the properties (statistical and other) of the received signal. 
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By considering the 'ratio' of the pdfs in Figs. HI El [H or [TT] relative to the Xg asymptotic 
distribution of the deep- field energy density, one obtains a transfer function for fs e (s e ) 
that captures the change in statistical properties, similar to the way in which Fresnel 
reflection or refraction coefficients relate incident to reflected or refracted fields. It also 
yields information on the likelihood and extent of signal distortion. 

• Measurement of constitutive EM properties of materials inside a mode-stirred reverberation 
chamber, including anisotropy, absorptivity, conductivity and/ or permittivity of materials: 
Several situations in which reverberant fields near a PEC surface are of practical relevance 
were already mentioned in [Arnaut and West, 2006]; here, we focus on media with finite 
o~/(ueo) and e only. Through the measurement of probability distributions of reflected 
or refracted fields at a fixed and sufficiently close distance from the boundary of a semi- 
infinite medium or, by extension, from a flat panel of sufficient large electrical thickness, 
constitutive parameter values can be deduced using Figs. HJ [9j [11] as a statistical 
inversion problem [Arnaut, 2006b]. Alternatively, the spatial functional fco^-dependence 
of the average electric energy density (particularly its normal component) can be compared 
with the results in Figs. [2] or [7] to deduce these parameters, circumventing the need for 
determining the pdf. 

7 Conclusion 

In this paper, we have been concerned with the influence of the presence of a deterministic 
semi- infinite isotropic medium on the local first-order statistical, i.e., distributional properties 
of an ideal random electromagnetic field for nondirectional (hemispherical) incidence. The EM 
boundary conditions cause statistical anisotropy of the vector field and a redistribution of the en- 
ergy density between its tangential and normal components, compared to the statistical isotropy 
and homogeneity of the incident field. This modifies its probability density function. For the 
reflected plus incident field, this modification changes in a damped oscillatory manner as a func- 
tion of distance from the interface; for refracted fields, the change is homogeneous within the 
dielectric because no standing waves occur in this region. It was found that, for the tangential 
and normal field components, the presence of a dielectric alters their variability (as expressed by 
their standard deviation) in an opposing manner, viz., for e/eo > 1 the amplitude fluctuations 
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of the refracted field become more compressed, whereas for the incident plus reflected field the 
fluctuations expand. 

In this paper, the analytical formulation in terms of numerically solved integrals was made 
possible by the simplicity of the geometry (single interface) and medium (isotropy). More 
complex configurations involving multiple, finite, angled or curved interfaces, e.g., for interior 
fields inside overmoded cavities or multiple-scattering problems, are unlikely to afford such an 
semi-analytical approach. In these cases, one may numerically calculate the statistics by using 
a Monte Carlo simulation for angular spectra of random plane waves (or other suitable random 
excitation) with specified input statistics. The numerically solved reflected and refracted fields 
can then be collated and evaluated at a single location to yield the output statistics for these 
fields (cf. Sec. V of [Arnaut and West, 2006]). 

Acknowledgment 

This work was sponsored in part by the 2003-2006 Electrical Programme of the UK Department 
of Trade and Industry National Measurement System Policy Unit (project no. E03E54). 

References 

[1] Arnaut, L. R. (2002), Compound exponential distributions for undermoded reverberation 
chambers, IEEE Trans. Electromagn. Compat., 44 (3), 442-457. 

[2] Arnaut, L. R. and West, P. D. (2006), Electromagnetic reverberation near a perfectly 
conducting boundary, IEEE Trans. Electromagn. Compat., ^S(2), 359-371. 

[3] Arnaut, L. R. (2006a), Spatial correlation functions of inhomogeneous random electromag- 
netic fields, Phys. Rev. E, 73(3), # 036604. 

[4] Arnaut, L. R. (2006b), Spatial correlation functions of random electromagnetic fields in the 
presence of a semi-infinite isotropic medium, Phys. Rev. E, 74(5), # 056610. 

[5] Booker, H. G. and Clemmow, P. C. (1950), The concept of an angular spectrum of plane 
waves, and its relation to that of polar diagram and aperture distribution, Proc. Instn. 
Electr. Engnrs., Pt. Ill, 97, 11-17. 



14 



[6] Bourret, R. C. (1960), Coherence properties of blackbody radiation, Nuovo Cimento, 
XVIII (2), 347-356. 

[7] Dunn, J. M. (1990), Local, high-frequency analysis of the fields in a mode-stirred chamber, 
IEEE Trans. Electromagn. Compat, 32(1), 53-58. 

[8] Eckhardt, B., Dorr, U., Kuhl, U., and Stockmann, H.-J. (1999), Correlations of electromag- 
netic fields in chaotic cavities, Europhys. Lett., ^#(2), 134-140. 

[9] Hill, D. A. (1998) Plane wave integral representation for fields in reverberation chambers, 
IEEE Trans. Electromagn. Compat., 40(3), 209-217. 

[10] Hill, D. A. and Ladbury, J. M. (2002), Spatial-correlation functions of fields and energy 
density in a reverberation chamber, IEEE Trans. Electromagn. Compat., 44(2), 95-101. 

[11] Mehta, C. L. and Wolf, E. (1964), Coherence properties of blackbody radiation. I. Correla- 
tion tensors of the classical field, Phys. Rev., 134 {5 A), A1143-A1149. 

[12] Sarfatt, J. (1963), Quantum-mechanical correlation theory of electromagnetic fields, Nuovo 
Cimento, XXVII (S), 1119-1129. 

[13] Whittaker, E. T. (1902), On the partial differential equations in mathematical physics, 
Math. Ann., LVII, 333-355, §§5.1, 5.2. 



15 




Figure 1: Coordinate system and local plane of incidence (<po = 0, 1^ = l y ) for single 
TE wave component reflected and refracted by a semi-infinite isotropic medium. 




Figure 2: Average energy densities for normal and tangential components of the 
incident plus reflected electric field at selected values of u/{oj€q), normalized with 
respect to the energy density of the incident vector field, (S eo ). Functions increasing 
at koz = 1 represent (S et ) / (S eo ); functions decreasing at k$z = 1 represent {S ez ) / '{S eo }. 
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Figure 3: Polarization coefficients (i = 1,2) for the electric energy density of the 
incident plus reflected field near a conducting medium, at selected values of a/(ueo). 
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Figure 4: Pdf of S e of the incident plus reflected field for selected values of a/(ueo) 
at koz = 7r/4. 



17 




S/(S) (dB) 



Figure 5: Pdf of S e of the incident plus reflected field for selected values of a/(toeo) 
at k z = 0.01. 
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Figure 6: Std of S e of the incident plus reflected total field: (a) as a function of k$z 
for selected values of cr/(weo); (b) as a function of a/{ioeo) at selected values of k$z. 
Dotted lines in figure (b) represent limit values for a PEC surface at the indicated 
values of koz. 
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Figure 7: Intensities of tangential and normal electric field components as a function 
of relative permittivity: (a) for the incident plus reflected field at koz = 20tt; (b) for 
the refracted field. 
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Figure 8: Polarization coefficients (i = 1,2) for a lossless dielectric medium: (a) 
for the incident plus reflected field, as a function of electrical distance from the 
boundary, at selected values of the permittivity; (b) for the incident plus reflected 
field (Rx; at arbitrarily large or small distance) and for the refracted field (Tx; at 
any distance) as a function of the relative^ permittivity. 
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Figure 9: Pdf of electric energy density S e of the refracted field for selected values 
of e r at arbitrary k$z. 



0.75 




Figure 10: Std of S e of the refracted field as a function of e r at arbitrary koz. 
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Figure 11: Pdf of S e of the incident plus reflected field for selected values of e r at 

k z = tt/128. 
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